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ABSTRACT 

We present the HERA CO-Line Extragalactic Survey (HERACLES), an atlas of CO emission from 
18 nearby galaxies that are also part of The H I Nearby Galaxy Survey (THINGS) and the Spitzer 
Infrared Nearby Galaxies Survey (SINGS). We used the HERA multi-pixel receiver on the IRAM 30-m 
telescope to map the CO J = 2 ^ 1 line over the full optical disk (defined by the isophotal radius 
^25) of each target, at 13" angular resolution and 2.6 km s~^ velocity resolution. Here we describe 
the observations and reduction of the data and show channel maps, azimuthally averaged profiles, 
integrated intensity maps, and peak intensity maps. The implied H2 masses range from 7 x 10^ to 
6 X 10^ Mq, with four low metallicity dwarf irregular galaxies yielding only upper limits. In the cases 
where CO is detected, the integrated H2-to-H I ratios range from 0.02 - 1.13 and H2-to-stellar mass 
ratios from 0.01 to 0.25. Exponential scale lengths of the CO emission for our targets are in the range 
0.8 - 3.2 kpc, or 0.2 ± 0.05 r25- The intensity- weighted mean velocity of CO matches that of H I very 
well, with a la scatter of only 6 km s~^. The CO J = 2 ^ 1/J=1^0 line ratio varies over a range 
similar to that found in the Milky Way and other nearby galaxies, ~ 0.6-1.0, with higher values found 
in the centers of galaxies. The typical line ratio, ^ 0.8, could be produced by optically thick gas with 
an excitation temperature of 10 K. 

Subject headings: galaxies: ISM — ISM: molecules — radio lines: galaxies 



1. INTRODUCTION 

Molecular hydrogen, H2 , is the phase of the interstel- 
lar medium (ISM) most closely related to star formation. 
Clouds of molecular gas are thought to be birthplaces 
of virtually all stars and H2 often dominates the mass 
budget of the interstellar medium (ISM) in the inner, 
most vigorously star-forming parts of spiral galaxies (< 
0.5 r25). Unfortunately, H2 lacks a dipole moment and 
typical temperatures in giant molecular clouds (GMCs) 
are too low to excite quadrupole or vibrational transi- 
tions. Therefore, indirect approaches are required to es- 
timate the distribution of molecular hydrogen. Although 
not without associated uncertainties, CO line emission 
remains the most straightforward and reliable tracer of 
H2 in galaxies. It is relatively bright and its ability to 
trace the bulk distribution of H2 has been confirmed via 
comparisons with gamma rays (e.g., iLebrun et aLlll983l: 
Strong et al.lll98^ and dust emission (e.g.. iDesert et al.l 
1988 1 iDame enLfl200l . 

Since the first detections of CO emission from Milky 
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Way molecular clouds (Wilson et al.' '1970) and other 
galaxies (Rickard ct al. 1975; Solomon & de Zafra 1975), 
a number of surveys have use d CO to character i ze the 
molecular content of galaxies. lYoung fc Scovilld ()199lD 
summarize the first two decades of such observations, 
during which CO was observed in ^ 100 galaxies. These 
data yielded a basic understanding of the radial distri- 
bution of CO in disk galaxies and the relationship be- 
tween H I (atomic hydrogen), H2, and star formation 
as a function of morphology. This number expanded 
to ~ 300 with the publication of The Five Colleges Ra- 
dio Astronomy Observatory (FCRAO) Extr agalactic CO 
Survey (henceforth, the "FCRAO survey," lYoung et al.l 
[T995[) . which remains the definitive survey of CO in the 
local volume. 

The FCRAO survey and most of its predecessors fo- 
cused on the CO J = 1 ^ transition, i.e. the tran- 
sition between t he first rotati onally excited state and 
the ground state. iBraine et al.l (fl993l) observed both CO 
J = 2 ^ 1 and J = 1 ^ emission from the centers of 
81 galaxies (finding a typical line ratio of 0.89 ± 0.06). 
These data support the idea that CO J = 2^1 emis- 
sion is a vi able tracer of H2 i n other galaxies: llsrael et al.l 
11984) and ISakamoto et all |l995) showed that the dis- 
tribution of CO J = 2 ^ 1 follows that of CO J = 1 ^ 
in the Milky Way. 

In the past decade, technical improvements — espe- 
cially the construction of several millimeter-wave inter- 
ferometers — have allowed higher resolution imaging 
of CO in galaxies. These have obtained maps of Lo- 
cal Group galaxies with spatial resolution matched to 
the sizes of individual GMCs (e.g. , Fukui c t al.l Il999l: 
Engargiola^sLall 120031 : iLerov et al.l l2006t iRosolowskvl 
2007) and observations of large samples of more distant 
galaxies with resolutions of only a few hundred parsecs 
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(e.g-. IWalter et al.ll2001l : ISakamoto et al.lll999D . Of par- 
ticular note, the Berkeley Illinois Maryland Association 
Survey of Nearby Galaxies (BIMA SONG; iHelfer et all 
[200I mapped CO emission at ~ 6" resolution in 44 
nearby spirals. The relatively wide extent of the BIMA 
SONG maps (typical full width ~ 190") make this sur- 
vey the first to systematically map the disks of normal 
spiral galaxies with resolution matched to the scales over 
which GMCs, the basic units of the molecular ISM, are 
believed to form. 

Over the last few years another technical improvement, 
heterodyne receiver arrays on single-dish telescopes ob- 
serving at 1-3 mm, has made it possible to efficiently map 
CO emission from large areas of the sky. This allows com- 
plete inventories of molecular gas in samples of nearby 
galaxies that were prohibitive with single-pixel receivers 
(the FCRAO survey , by contrast, only sam pled the major 
axes of its targets). ISchuster et al.l ()2007D demonstrated 
this cap ability by using the Heterodyne Receiver Array 
(HERA. ISchuster et al.ll2004f l on the IRAM 30-m to map 
CO J = 2 — » 1 em ission from the whole disk of M 51. 
iKuno et al.l (|2007t ) presented CO J = 1 ^ maps of 40 
nearby spirals, largely carried out using the BEARS fo- 
cal plane array on the Nobeyama 45-m. In both cases, 
the large diameter of the telescopes (and the fact that 
HERA observes at 230 GHz) ensure a resolution of 11- 
15", adequate to resolve spiral structure, bars, rings, and 
large star-forming complexes. 

Here we present the HERA CO-Line Extragalactic 
Survey (HERACLES). HERACLES is a new survey of 
CO emission in nearby galaxies carried out with the 
HERA receiver array on the IRAM 30-m telescope^^. 
The main motivation for this new survey is to quantify 
the relationship between atomic gas, molecular gas, and 
star formation in a significant sample of galaxies. To 
meet this goal, HERACLES differs from previous CO 
surveys in two main ways: the area surveyed and the 
sample. 

In contrast to the FCRAO survey or BIMA SONG, 
we map the whole optical disk of each galaxy with full 
spatial sampling and good sensitivity. Despite the large 
area surveyed, we still achieve a relatively good resolu- 
tion of 13", ~ 500 pc at the median distance of our sam- 
ple. Because we use a single dish telescope rather than 
an interferometer, we are sensitive to extended structure 
that may be missed by the latter (which do not recover 
extended structure on the sky by design). 

The HERACLES sample ove rlaps that of The 
HI Nearby Galaxy Survey ([Walter et all l2008f) 
and the Spitzer Inf rared Nearby Galaxies Survey 
(jKennicutt et al.ll2003f l. ensuring an immediately acces- 
sible multiwavelength database spanning from radio to 
UV. These data offer windows on the diffuse ISM, dust, 
embedded star formation, photodissociation regions, 
and both young and old stars. 

This paper describes the HERACLES observing (© 
and data reduction (Sj3]) strategies. We show the distri- 
bution of CO emission and compare HERACLES 
results to previous observations of the H I 21-cm fi j5.3p 
and the CO J = 1 ^ f !j5.4p lines in the same galaxies. 
Two papers present the first round of in-depth scientific 

" IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the ION (Spain) 



TABLE 1 
HERACLES Sample 



Galaxy 


Dist. 


Incl. 


P.A. 


b 

f25 




(Mpc) 


(°) 


(°) 


(') 


NGC fi2S 


7.3 


7 


20 


4.9 




9.2 


66 


287 


5.3 


HolmbGrg II 




41 


1 77 


0.0 


NOP 9J^41 


14.1 


74 


154 


5.3 


NGC 2Qn'^ 


8.9 


65 


204 


5.9 


T— Tr^l Tvi r^£>i* rr T 


3.8 


12 


50 


1.6 


NGC 2976 


3.6 


65 


335 


3.6 


NGC 3184 


11.1 


16 


179 


3.7 


NGC 3198 


13.8 


72 


215 


3.2 


IC 2574 


4.0 


53 


56 


6.4 


NGC 3351 


10.1 


41 


192 


3.6 


NGC 3521 


10.7 


73 


340 


4.2 


NGC 4214 


2.9 


44 


65 


3.4 


NGC 4736 


4.7 


41 


296 


3.9 


DDO 154 


4.3 


66 


230 


1.0 


NGC 5055 


10.1 


59 


102 


5.9 


NGC 6946 


5.9 


33 


243 


5.7 


NGC 7331 


14.7 


76 


168 


4.6 



Position angle of major axis, measured 
north through east.*^ Radius of the B-band 
25 mag/arcsec^ isophote. 
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Fig. 1. — Our observing strategy illustrated for the inclined 
spiral galaxy NGC 3521. The gray scale ima ge shows the THINGS 
H I column density map llWalter et al.l l2008) and the oval indicates 
r25, our target area. Arrows indicate the length and direction of 
individual scan legs, which form back— and— forth pairs that together 
yield a fully-sampled map. To cover the galaxy we use 7 back- 
and-forth scans along the short edge and 3 along the long edge 
(indicated by the entry 7 X 3 in column 3 of Table [2J. The beam 
print of HERA appears at the reference position, which was chosen 
to be free of gas but near the galaxy. Concentric rectangles show 
the areas of full (inner) and partial (outer) sensitivity. 

analysis. In lBigiel et aH ()2008l ). we combine HERACLES 
with Hi, IR, and UV data to measure the relationship 
between the surface densities of H I , H2, and star for - 
mation rate (the "Schmidt Law"). In lLerov et al.l ()2008f ) 
we use the same data to measure the star formation per 
unit gas as a function of environment, comparing these 
measurements to proposed theories. 

2. OBSERVATIONS 
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TABLE 2 
Observing Dates and Orientation 





Scan P.A.'' 




DflytrGs Observed 


Spectra' 




\ ) 






(iU ) 


NGC 628 





9x9 


10, 12, 13, 14 Jan 2007 


105 


NGC 925 


287 


9x4 


23^*, 30 Oct 2007; 12°, 13 Jan 2008; 25, 26 Mar 2008 


55 


Holmberg II 





5x5 


21, 22 Feb 2008 


28 


NGC 2841 


153 


6x2 


26, 27, 30 Oct 2007 


31 


NGC 2903 


204 


10 X 4 


27'' Nov 2007; 13, 16, 17, 21 Feb 2008*= 


87 


Holmberg I 





3x3 


30 Oct 2007; 1, 3 Nov 2007 


18 


NGC 2976 


335 


6x3 


27, 28, 31 Oct 2007; 3 Nov 2007 


44 


NGC 3184 





6x6 


31 Oct 2007; 1, 4, 5 Nov 2007 


66 


NGC 3198 


215 


6x2 


26, 28 Oct 2007; 2 Nov 2007 


31 


IC 2574 


55 


11 x 6 


17, 21, 22, 25, 27 Feb 2008; 25, 26 Mar 2008 


125 


NGC 3351 





6x6 


31 Oct 2007; 1, 2 Nov 2007; 17 Jan 2008 


70 


NGC 3521 


340 


7x3 


2, 3, 4 Nov 2007 


50 


NGC 4214 


65 


4x4 


5 Nov 2007; 10, 11, 12 Jan 2008 


38 


NGC 4736 





8x4 


19, 20 Jan 2007; 21, 24<^, 25, 26 Feb 2007 


75 


DDO 154 


230 


4x2 


13, 16 Feb 2008 


16 


NGC 5055 





13 x 5 


13°, 14, 15, 16 Jan 2007 


95 


NGC 6946 





11 X 11 


13, 14, 15, 17, 18, 19 Jan 2007 


207 


NGC 7331 


168 


8x3 


23^*, 24^^, 26, 27 Oct 2007; 10, 11, 12 Jan 2008 


63 



^ Orientation of major axis scans measured from north through east. See Figure[T]'' Number of fully- 
sampled (i.e., back and forth) scans along short x long edge. Width of one scan is ~ 68".° Number 
of spectra in the final reduction. Each represents an 0.5 second integration by one polarization on 
one receiver. One hour integrating on source yields ~ 13 x 10** spectra.'^ For this day, we adjusted 
the rejection th resho ld to remove more than the 10% of spectra with the highest RMS about the 
baseline fit (see >i3.3|l .° For this day, only one polarization was usable. 



The HERACLES sample consists of galaxies that are 
targets of THINGS, far enough north to be easily ob- 
served by the 30-m, and not prohibitively large (< 12 x 
12'). These 18 spiral and irregular galaxies are listed in 
Table [T] We draw the inc lination, position an gle, and 
distances in this table from I Walter et ahl (|2008f) and the 
isophotal radius, r25, from the LEDA database. 

From January 2006 through March 2008, we used the 
HERA on the IRAM 30-m telescope to map CO J = 
2-^1 emission from these galaxies. HERA is a multi- 
pixel receiver that simultaneously observes 9 positions on 
the sky at two orthogonal linear polarizations. 

We tuned HERA to observe the CO J = 2 ^ 1 transi- 
tion near 230 GHz and attached the HERA receivers to 
the Wideband Line Multiple Autocorrelator (WILMA) 
backend. WILMA consists of 18 units, each with 2 MHz 
channel width and 930 MHz bandwidth, yielding a ve- 
locity resolution of 2.6 km s^^ and velocity coverage of 
1200 km s-i at the wavelength of the CO J = 2 ^ 1 
line. 

Table [2] (column 4) gives the dates when individual 
galaxies were observed. We generally observed only un- 
der good winter conditions, meaning zenith opacities 
T < 0.3 at 225 GHz. During most runs, the opacity 
was better than this, often < 0.1. System temperatures 
varied with pixel, receiver, and semester but typical val- 
ues were 275-350 K for the 9 receiver pixels of the first 
polarization (HERAl) and 350-450 K for the 9 receiver 
pixels of the second polarization (HERA2). On average, 
data from one pixel (not always the same one) could not 
be used due to unreliable baselines. During a few days, 
indicated in Tabic [2l only one polarization was usable. 

We used on-the-fly mapping mode and scanned across 
the galaxy at 8" sec~^, writing out a spectrum every 0.5 
seconds, i.e., integrating over 4" of scanning. Approxi- 
mately every 2 minutes, we observed a reference position 



near the galaxy for 10-15 seconds. This was chosen us- 
ing the THINGS H I maps to pick a position free of gas 
that lies outside but near the optical disk of the galaxy. 
Every 10-15 minutes we calibrated the intensity scale of 
the data using sky emission from the reference position 
and the backend counts from loads at ambient and cold 
(liquid nitrogen) temperatures. 

In most cases the target area was the optical disk of the 
galaxy, defined by the 25**^ magnitude i?-band isophote 
('"25), though in several galaxies we extended the map to 
probe obvious H I peaks or filaments outside the optical 
disk. 

To construct a single map of a galaxy, we observed the 
target area using a series of parallel scans. These "long 
edge" scans were usually aligned with the major axis 
of the galaxy and were spaced to cover the whole target 
area. We then immediately made a series of "short edge" 
scans covering the same area, but oriented perpendicular 
to the original long edge scans. The goal of this cross- 
hatched observing strategy was to minimize artifacts in 
the final data by ensuring that each area of the sky was 
observed by many different receivers. Figure [1] illustrates 
this approach for NGC 3521, an inclined Sb galaxy. The 
beam-print of the array is shown at the reference po- 
sition and arrows over the body of the galaxy indicate 
individual scans. Here there are 3 "long edge" scans and 
7 "short edge" scans. The second and third columns in 
Table[2]give the orientation of the long edge — measured 
north through east — and the number of scans along the 
short and long edges (7x3 for NGC 3521). 

In order to fully sample the 11" (FWHM) beam of the 
30-m, we used the derotator to rotate the beam pattern 
of HERA by 9.5° relative to the direction of scanning^^. 

For a full description and illustration of this observing mode 
see the HERA User Manual (Version 2) by Schuster et al. (2006), 
available online at ,http : //www. iram. f r/IRAMES| 
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This yields a pixel spacing of ~ 4" on the sky but leaves 
a gap between the HERA pixels, which are separated 
by 28". We filled this gap by offsetting the next scan 
by ff.9" perpendicular to the scan direction and then 
repeating the original scan leg in reverse. This strategy 
yields a fully sampled map with spectra spaced by 4" 
both along and perpendicular to the scan direction. 

We repeated the full set of long and short edge scans 
15-20 times for each galaxy. This yielded the equivalent 
of 2-3 minutes of integration per independent beam. The 
fifth column in Table [2] lists the number of individual 
spectra used in the final data cube for each galaxy (in 
units of 10''). Each spectra corresponds to a 0.5 second 
integration with one receiver, so that the full 18-element 
array produces 36 spectra each second of on-source inte- 
gration and an hour on source produces 13 x lO'* spectra. 
Note that these numbers come after rejection of high- 
RMS spectra and removal of bad pixels and that 
this time estimate takes no account of reference observa- 
tions, calibrations, slewing, tuning, or other overheads. 

3. REDUCTION 

3.1. Basic Calibration and Reference Subtraction 

We carried out the basic data reduction using the Mul- 
tichannel Imaging and Calibration Software for Receiver 
Arrays (MIRA)^'^, which is part of the Grenoble Image 
and Line Data Analysis Software (GILDAS) package^"*. 
This consisted of combining each scan with the nearest 
reference measurement and calibration observation via 



Ta — Teal X 



On -Off 
Hot - Off 



(1) 



Here On refers to the backend count rate from the on- 
source measurement. Off to the count rate from the ref- 
erence (empty sky) measurements before and after the 
on-source measurement. Hot to the count rate from the 
ambient load during the calibration observation, and Tcai 
to a calibration factor determined from the calibration 
observation^^. The resulting antenna temperature, T^, 
is defined to be the brightness temperature of a source 
filling the entire 2tt steradians in front of the telescope 
and outside the atmosphere. 

The data, in units of T^, were written out for further 
reduction using the Continuum and Line Analysis Single- 
dish Software (CLASS) package, which is also part of 
GILDAS. 

3.2. Fitting Baselines Based on H I Velocity 

After this basic reduction, total power variations in 
both the receiver and the atmosphere and nonlinearities 
in the backend and receiver still made it necessary to 
fit baselines to individual spectra before combining them 
into data cubes. Because the data total > 10'' spectra, 
this needed to be done in an automated manner. 

The CO line typically covers a small but variable part 
of the bandpass, so fitting a linear baseline to a single 
wide window did not yield satisfactory results. Instead, 
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to the IRAM 30-m telescope, see "Calibration of Spectral Lino Data 
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we used the THINGS H I data cubes to define the re- 
gion of the spectrum likely to contain the CO line and to 
fit a linear baseline to a restricted part of the spectrum. 
This yielded good results and was straightforward to au- 
tomate. The underlying assumption is that H I and H2 
(traced by CO) are reasonably well-mixed, so that the 
mean velocity of CO emission is similar to the mean H I 
velocity. This is a reasonable expectation: the two lines 
both trace dissipative gas moving in the same potential 
well and H2 is believed to form out of H I. In ^5.'S\ we 
use our data to verify this assumption, indeed finding an 
excellent match between the mean velocities of H I and 
CO. 

The exact approach, illustrated in Figure [21 was: 

1. We defined a window over which the CO line was 
likely to appear. This window was centered on the 
mean H I velocity and 60-300 km s~^ wide, de- 
pending on the galaxy. For many galaxies, we de- 
fined two windows: one used near the center of the 
galaxy and another used over the rest of the disk 
with the fitting region varying smoothly between 
the two regimes. We selected the window widths 
for each galaxy based on a preliminary reduction. 

2. Adjacent to this window, we defined two regions of 
the spectrum that we used to fit a linear baseline. 
These had the same width as the central window, 
which was not included in the fit. The fit is sub- 
tracted from the whole spectrum. 

3. We blank all data outside the fitting windows, so 
that only the fitting windows and the central region 
— the likely location of signal — are left. 

In this way, we fit baselines to all spectra in each 
galaxy. We reduced each data cube three times. First, 
we made a crude reduction fitting baselines to a fairly 
wide area around the H I line. Based on the results, we 
refined the fitting region (white areas in Figure [2l step 
1 above), defining a central window that was as narrow 
as possible while still including all of the CO emission 
seen in the preliminary reduction. With the third reduc- 
tion, we iterated the process, using the second reduction 
to further refine our fitting regions and more carefully 
identifying and removing pathological data (i 33.3|) . 

3.3. Rejecting Remaining Pathological Spectra 

A small fraction of spectra are still not well-fit by our 
approach. If included in the final data cube, these in- 
troduce artifacts and obscure signal. To remove them in 
a straightforward, systematic way, we discard the 10% 
of the data with the highest RMS residuals about the 
baseline fit from each day of observing (note that these 
residuals are determined from the fitting region, which 
is chosen to be signal free). Removing 10% of the data 
causes a negligible decrease in signal-to-noise but remov- 
ing the pathological spectra yields a noticeable increase 
in the quality of the data cubes. In a few cases where a 
htife ar population of high RMS spectra survive the baseline 
Hlting, we adjust the rejection threshold to reject more 
data (these observations are indicated in Table [2|) . Typ- 
ical values for the 90*'^ percentile are RMS (TJ) ~ 0.35- 
0.7 K for 2.6 km s~^ wide channels. 
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Fig. 2. — Our approach to baseline fitting, illustrated for a single (0.5 second) spectrum of NGC 6946. We use the mean H I velocity to 
predict the likely velocity of CO emission. From this, we define three regions: one near the line that is likely to contain signal and two that 
are close to the line but displaced enough not to contain signal. We use the data in these two adjacent regions to fit a first-order baseline 
(while avoiding the central region). We subtract this baseline from the whole spectrum and discard data outside the fitting regions. 



3.4. Constructing a Data Cube 

We combine the reduced spectra into a table and use 
the CLASS gridding routine xy_map to construct a data 
cube with 2" pixel size and 2.6 km s""'^ channel width. 
This routine convolves the irregularly gridded OTF data 
with a Gaussian kernel with a full with width of ^ 1/3 
the FWHM beam size, yielding a final angular resolution 
approximately 11" and noise correlated on scales of 4". 

Because the noise is correlated on scales that are small 
compared to the beam (response to astronomical signal), 
smoothing the data slightly offers a significant gain in 
sensitivity while min imally degrading the beam size (e.g., 
iMangum et al.|[2007[ ). In the rest of this paper, we show 
data that have been convolved with a 7" (FWHM) Gaus- 
sian, yielding a final resolution of 13". 

We converted the units of the data cube from antenna 
temperature, TJ, to main beam temperature, Tmb- Main 



beam temperature is the temperature of a source filling 
only the main beam of the telescope (recall that refers 
to a source filling the full forward 2tt of the sky). To con- 
vert from to Tmb, we replaced the forward efficiency, 
Fcs ~ 0.91, with the main beam efficiency, BcB — 0.52, 
yielding Tmb = 0.91/0.52 = 1.75 (efficiencies were 
adopted from the 30-m online documentation). For the 
remainder of the paper, we work in units of Tmb- 

3.5. Uncertainties 

Table[3]lists the RMS noise for each galaxy (column 2). 
We estimate this from signal-free parts of the cube, which 
has 13" angular resolution and 2.6 km s~^ channel width. 
Typical values are in the range 20 - 25 mK. Mostly this 
noise averages in the expected manner, but on very large 
scales our method of calibration introduces an additional 
consideration. Because we subtract the same reference 
measurements from all spectra along an individual scan 
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leg (Section[2]), low level correlated noise is present in the 
maps. As a result, the RMS noise of spectra derived from 
integrating entire data cubes is 2-3 times higher than 
expected from averaging independent, normal noise; this 
ratio is consistent with the ratio of observing time spent 
on the reference to that on source. 

We calculated how the integrated intensity of regions 
with very high (> 50) signal-to-noise ratio (SNR) var- 
ied when measured with different polarizations and on 
different days. We identified such regions in NGC 2903, 
3351, 3521, 5055, and 6946 and then reduced the data 
separately for each day and polarization. The RMS day- 
to-day scatter in the integrated intensity from these high 
SNR regions is ^ 20%. This is an approximate, but 
by no means rigorous, measure of the uncertainty in the 
calibration of the telescope (e.g., it includes the effect of 
pointing errors but not the uncertainty in the assumed 
efficiencies of the telescope). 

We tested the uncertainty associated with our method 
of reduction and baseline fitting by reducing all of the 
data for one galaxy independently. We flagged bad data 
and identified baseline fitting regions by eye. Over re- 
gions of significant emission (/co > 2(t in 3 successive 
channels at 30" resolution) the two reductions agree very 
well. The best-fit line relating our automated reduction 
to the "by-hand" one in these regions has a slope of 1.01 
and an intercept of 0.06 K km s~^. Along lines of sight 
with strong CO emission (/co > 2.6 K km s~^), the ratio 
of the two reductions scatters by 10% (RMS). This check 
offers an estimate of the uncertainty associated with how 
we carry out the data reduction and verifies that our 
automated approach (which is less subjective and much 
simpler to apply) agrees well with a "by-hand" reduction. 

A final important uncertainty is related to our observ- 
ing strategy. In some cases, pixel-to-pixel variations in 
gain, system temperature, and bandpass shape leave the 
imprint of our cross-hatch observing strategy in the data. 
This may be seen in the channel and integrated intensity 
maps present (21) as striping along the two scan direc- 
tions (i.e., it is possible to see the scanning path of indi- 
vidual pixels in the noise) . The magnitude of this striping 
is usually comparable to or below that of the statistical 
noise and mostly these artifacts do not obscure or mimic 
signal. The exceptions are four highly inclined spirals — 
NGC 2841, 2903, 3521, and 7331 — in which significant 
striping is visible in individual channel maps. The likely 
cause is a breakdown in our baseline fitting: because the 
CO line is often wide in these galaxies, we are forced to 
use relatively broad fitting windows. 

4. THE DISTRIBUTION OF CO J = 2 ^ 1 EMISSION 

Our basic result is the distribution of CO J = 2 ^ 1 
emission in our 18 targets. In 14 galaxies, extended emis- 
sion is clearly detected. Here, we show that emission 
using channel maps, peak intensity maps, integrated in- 
tensity maps, and radial profiles. We make several basic 
measurements, summarized in Table [U including the in- 
tegrated fiux, luminosity, and exponential scale length. 

4.1. Channel, Integrated, and Peak Intensity Maps 

Figures [3l - (TBI show channel maps (i.e., intensity inte- 
grated over a succession of velocity ranges) of each clearly 
detected galaxy. These span approximately the velocity 
range over which we find CO emission and show the area 



of uniform sensitivity (i.e., the inner rectangle in Fig- 
ure [T]) . The gray scale and final channel width vary from 
galaxy to galaxy and are indicated in each caption. Data 
at velocities outside the baseline fitting region f !j3.2p are 
blanked in the reduction and appear as zeroes (smooth 
gray regions) in the channel maps. 

The left column in Figures [T7] - [301 shows integrated 
intensity for the same 14 galaxies. A bar next to each fig- 
ure indicates the color stretch, which varies from galaxy 
to galaxy. Contours in the integrated intensity maps 
run from 1 K km s~^ (gray), increasing by factors of 
2 each step. The first black contour, /co = 2 K km s~^, 
corresponds to Eh2 ^11 M0 pc~^ in a face-on spiral 
galaxy (Equation!?]) . This approximately indicates where 
the ISM is dominated by H2. The first white contour, 
/co = 32 K km s^^, corresponds to ^ 176 M© pc^^, 
roughl y the surface densit y of an individual Galactic 
CMC (|Solomon et al.l[l98l . Re gions with such high sur- 
face densities over a large area are rare but not unknown 
in our survey (e.g., the bar of NGC 2903 and the centers 
of NGC 2903, 3351, 4736, 5055, 6946). 

We created the integrated intensity maps by summing 
along the velocity axis over regions that contain signal. 
To identify these regions, we first smooth the data to 
30" angular resolution and bin to reach a velocity reso- 
lution equal to half the width of the channel maps (this 
varies with galaxy, see Figures |3l - fTB]) . All regions with 
/co > 3(7 in two consecutive velocity channels in these 
low resolution cubes were labeled "signal," and we inte- 
grated the original cube over these regions to yield the 
maps in Figures [T7] - [301 This procedure has a low false 
positive rate (< 1 expected per cube) so that most emis- 
sion seen here is real with high confidence, though in a 
few cases (NGC 2841, 2903, 7331) artifacts from baseline 
fitting persist. Because the condition to identify signal 
is fairly restrictive these maps are not ideal indicators of 
low brightness, extended emission. 

The peak intensity maps in the right column of Figures 
[T7l-[30lgive a clearer view of the area mapped, extent of 
low brightness emission, and morphology. These are cre- 
ated by measuring the peak brightness along each line of 
sight at 2.6 km s~^ velocity resolution. This procedure 
suppresses the worst artifacts in our data, which tend to 
persist at a low level over many consecutive channels. It 
also lowers the contrast between bright and faint regions, 
which arise partially from differences in the line width. 
As a result, several interesting low-lying features are vis- 
ible in these maps: e.g., arms extending out from the 
molecular rings in NGC 3521 and NGC 7331; faint fila- 
mentary structure in the disk of NGC 4736; and a faint, 
previously unidentified molecular complex in the south- 
west of the dwarf starburst NGC 4214 fc.f.. [Walter et all 
[2OOII) . 

Peak intensity maps of inclined systems tend to show 
radial structure on small scales. In order for emission 
to add constructively in such a map, it must be at the 
same velocity. When a galaxy with a rapidly changing 
velocity field is observed with finite angular resolution 
and then collapsed to a peak intensity map, signal will 
tend to be smoothed along isovelocity contours, leaving 
them faintly visible in the final image. 

The channel, integrated, and peak intensity maps show 
extended molecular structures covering the inner part 
of the disk in most of our galaxies. The detailed mor- 
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Fig. 3. — Channel maps for NGC 628 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 5.2 km s^^ wide 
and the RMS intensity in one channel is 19 mK. The central LSR velocity for each channel (in km s~^) appears in the top left corner of 
the map. 
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Fig. 4. — Channel maps for NGC 925 at 13" resolution. The grayscale runs from -0.025 to 0.10 K. Each channel is 13 km s~^ wide and 
the RMS intensity in one channel i s 12 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( i|3.2l l. 

phologies are fairly varied. Spiral structure is partic- 
ularly evident in the Sc galaxies NGC 628, NGC3184, 
and NGC 6946. NGC 2903 is dominated by a central, 
bright bar. The Sb galaxies NGC 3351 and NGC 7331 
both show well-defined molecular rings and similar struc- 
tures are suggested by the channel maps for NGC 2841 
and NGC 3198. Molecular gas completely covers the 
inner parts of NGC 3521, NGC 4736, and NGC 5055. 
Over this area, NGC 3521 and NGC 4736 show ring- 
like structure and outside this central region, NGC 3521, 
NGC 4736, NGC 5055, and NGC 7331 show spiral struc- 
ture. The late-type galaxies NGC 925 and NGC 4214 
both show only faint CO emission, which is common for 
late-type galaxies. The other low mass galaxy that we 
detect, NGC 2976, defies this trend, showing bright CO 
dominated by two large star-forming complexes. 



4.2. Radial Profiles 



Despite different morphologies, our targets show a 
common behavior on large scales: azimuthally averaged 
CO intensity declines steadily as a function of radius in 
most galaxies. This is clear from Figure [311 in which 
we plot integrated intensity as a function of radius for 
each detected galaxy. To make these profiles, we collapse 
the data cubes without masking into integrated intensity 
maps and then average these over a series of concentric 
tilted rings. These rings assume the position angle and 
inclinations in Tabled] and are each 10" in width, about 
1 resolution element. Note that these profiles do not in- 
clude any correction for inclination (i.e., they correspond 
to the observed values). Error bars show the uncertainty 
in the mean integrated intensity, derived by measuring 
the RMS scatter within the tilted ring and dividing by 
^/N, where N is the number of independent data points 
in the ring. Measurements with > 3cr significance appear 
as black circles, those with less as gray triangles. 



Leroy et al. 




iff H" Ka"^w" kT 



Fig. 5. — Channel maps for NGC 2841 at 13" resolution. The grayscale runs from -0.025 to 0.075 K. Each channel is 65 km wide 
and the RMS intensity in one channel is 13 m K. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( ^3.2| I. 



TABLE 3 



Basic Results 


Galaxy 


CT rms 




Leo 




dmax 


A/h2/Mhi 


A/H2/M, 


{AfH2 +Mhi)/M, 




(mK) 


(10^ K km s ^ arcsec^) 


(10^ K km s-i pc^) 


(kpc) 


{r25) 








NGC 628 


22 


1.8 


23 


2.4 


0.72 


0.24 


0.10 


0.52 


NGC 925 


21 


0.25 


4.9 


3.2 


0.17 


0.04 


0.04 


0.90 


Holmberg II 


36 


< O.IO^' 


< 0.33=^ 










> 3.6 


NGC 2841 


44 


0.61 


28 




0.65 


0.13 


0.03 


0.22 


NGC 2903 


23 


2.9 


53 


1.6 


0.47 


0.50 


0.25 


0.77 


Holmberg I 


24 


< 0.34^^ 


< 0.94=^ 










> 7.2 


NGC 2976 


21 


0.40 


1.2 


1.2 


0.53 


0.36 


0.05 


0.18 


NGC 3184 


20 


1.0 


30 


2.9 


0.55 


0.39 


0.09 


0.33 


NGC 3198 


17 


0.25 


11 


2.7 


0.43 


0.04 


0.04 


1.03 


IG 2574 


33 


< 0.74=^ 


< 2.8=^ 










> 3.8 


NGC 3351 


20 


0.78 


19 


2.5 


0.41 


0.63 


0.05 


0.12 


NGC 3521 


22 


2.8 


76 


2.2 


0.95 


0.38 


0.08 


0.29 


NGC 4214 


18 


0.09 


0.17 




0.33 


0.02 


0.01 


0.80 


NGC 4736 


23 


2.2 


11 


0.8 


1.33 


1.13 


0.03 


0.06 


DDO 154 


18 


< 0.05=^ 


< 0.21=^ 










> 34.8 


NGC 5055 


24 


4.1 


98 


3.1 


0.67 


0.43 


0.08 


0.26 


NGC 6946 


25 


10.8 


88 


1.9 


0.86 


0.86 


0.16 


0.35 


NGC 7331 


20 


2.2 


113 


3.1 


0.61 


0.50 


0.07 


0.21 



Note. — Column 1: galaxy; Column 2: RMS intensity n a 2.6 km s ^-wide channel map with 13" (FWHM) angular 
resolution; Column 3: integrated flux of CO emission (see Equation [2] for conversion to Jy km s~^); Column 4-' integrated 
luminosity of CO emission (see Equation |3] for conversion to H2 mass); Column 5: exponential scale length of CO emission; 
Column 6: radius of most extended high SNR detection ( <^4.2|l : Column 8: ratio of H2 to H I mass; Column 9: ratio of H2 
to stellar mass; Column 10: ratio of total gas (H2 + H I) to stellar mass. 
^ 5a upper limits. 
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Fig. 6. — Channel maps for NGC 2903 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 28.6 km s~^ wide 
and the RMS intensity in one channel is 10 m K. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( i|3.2l l. 
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Fig. 7. — Channel maps for NGC 2976 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 20.8 km s"^ wide 
and the RMS intensity in one chan nel is 9 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( ^3.2| |. 
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Fig. 8. — Channel maps for NGC 3184 at 13" resolution. The grayscale runs from -0.025 to 0.10 K. Each channel is 15.6 km wide 
and the RMS intensity in one channel is 11 m K. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( i|3.2l l. 



Thick gray lines show exponential fits to the profiles, 
carried out over the region bracketed by the arrows at 
the top of each plot. The corresponding scale lengths, 
given in Table [31 span from 0.8 to 3.2 kpc. In most 
galaxies, an exponential decline captures the large-scale 
behavior well, though we cannot reliably fit NGC 2841 or 
NGC 42 14. Studying BIMA SONG, which overlaps our 



sample, Regan et al 



g mMA 
T (I200I 



found the same variations 



in detailed morphology averaging into a common expo- 
nential decline. They discuss this topic in more detail. 



Both iRegan et al] (|200ir ) and iHelfer et al.) (|2003D noted 
deviations from these fits in galaxy centers, which we 
also observe in the form of depressions (e.g., NGC 3521, 
NGC 7331) and excesses (e.g., NGC 3351, NGC 6946). 

The exponential decline results from the combination 
of a decreasing filling fraction of CO emission and de- 
creased intensity along individual lines of sight. To illus- 
trate the relative contributions of these two effects. Fig- 
ure [3T] also shows the maximum integrated intensity in- 
side each ring. Filled diamonds indicate rings where this 
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Fig. 9. — Channel maps for NGC 3198 at 13" resolution. The grayscale runs from -0.025 to 0.10 K 
and the RMS intensity in one chan nel is 4 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( ^3.2| I. 
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Fig. 10. — Channel maps for NGC 3351 at 13" resolution. The grayscale runs from -0.025 to 0.08 K. Each channel is 33.8 km s"! wide 
and the RMS intensity in one chan nel is 7 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( ^3.2| I. 



maximum value is > Aa and significantly (> 2a) larger 
than the absolute value of the minimum integrated inten- 
sity (roughly accounting for artifacts). Open diamonds 
show where the maximum intensity is below this value 
and we are correspondingly uncertain that the measure- 
ment is not biased by noise and artifacts. Comparing the 
two profiles, we observe a fairly close correspondence in 
the central parts of galaxies and then a shallow decline 
in maximum intensity compared to average intensity at 
larger radii. The difference corresponds to the decreasing 
area subtended by bright CO emission. 

Although H I typically dominates the ISM outside 
~ 0.5 r25, we detect individual regions in the outskirts 
of several galaxies at high significance. In Figure [?T] and 
Table [3l we report the largest radius at which we detect 
a "high signal spectrum." We define this as a region with 
three consecutive velocity channels at > 3tT significance 
over 10 spatially contiguous pixels. This conservative 
definition generates no false positives in the whole sur- 
vey when applied to the negative portions of the data. 
Figure [21] clearly shows that more extended, but lower 
significance, signal is present in many maps. 

For the most part, these detections are simply associ- 
ated with the outer part of the star forming disk and can 



be clearly identified with structures seen in H I, UV, or 
IR maps that extend continuously from the inner to outer 
galaxy. This is the case for the detections at ~ 0.7-1.0 r25 
in NGC 628, NGC 3521, NGC 5055, and NGC 6946 (the 
outer region in the latter corresponds to po sitions PIO, 
11, and 12 observed bv iBraine et all 120071 ). The no- 
table exception is the bright CO emission at 1.33 in 
NGC 4736. This relatively narrow-line width feature is 
clearly visible in the channel map centered at 405 km 
in Figure [42l The emission is ^ 5' (projected distance of 
~ 7 kpc) from the center of the galaxy and well outside 
the main region o f active star formation (e.g., see maps in 
iLerov et al.ll2008h . This feature matches the velocity and 
position of an H I filament outside the main body of the 
galaxy; it also shows UV and mid-IR (24/xm) emission, 
indications of recent massive star formation. 



4.3. Flux, Luminosity, and H2 Mass 

Columns 3 and 4 of Table [3] report the integrated CO 
flux and luminosity for each target, measured by sum- 
ming the same integrated intensity maps used to make 
the radial profiles. In all cases, the statistical uncertainty 
from Gaussian noise alone is small, meaning that the 
choice of area to integrate over, artifacts, and calibration 
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Fig. 11. — Channel maps for NGC 3521 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 41.6 km s^^ wide 
and the RMS intensity in one channel is 9 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( i\'i.'2l . 
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Fig. 12. — Channel maps for NGC 4214 at 13" resolution. The grayscale runs from -0.025 to 0.1 K. Each channel is 13.0 km s"'^ wide 
and the RMS intensity in one channel is 10 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( !)3.2|l . 
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Fig. 13. — Channel maps for NGC 4736 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 18.2 km wide 
and the RMS intensity in one channel is 12 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( i|3.2l l. 
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Fig. 14. — Channel maps for NGC 5055 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 28.6 km s^^ wide 
and the RMS intensity in one channel is 10 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( i|3.2l l. 
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Fig. 15. — Channel maps for NGC 6946 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 18.2 km s~^ wide 
and the RMS intensity in one channel is 13 m K. The smooth gray area indicates data at velocities outside the baseline fitting window, 
which are blanked during the reduction ( i|3.2l l. 
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Fig. 16. — Channel maps for NGC 7331 at 13" resolution. The grayscale runs from -0.025 to 0.15 K. Each channel is 52.0 km s~^ wide 
and the RMS intensity in one chan nel is 8 mK. The smooth gray area indicates data at velocities outside the baseline fitting window, which 
are blanked during the reduction ( ^3.2| I. 




Fig. 17. — {left) Integrated intensity CO J= 2 — > 1 intensity along each line of sight, summing over regions where Iqo > 2cr in 
consecutive channel maps (see text). Contours begin at 1 K km and increase by a factor of two each step. The first black contour 
indicates 2 K km s~^ and the first white contour shows 32 K km s~^, which correspond approximately to surface densities of 11 and 
176 Mq pc~^. In both figures a dashed ellipse shows r25, the 25*'' magnitude isophote, and a solid black circle indicates the 13" angular 
resolution of the data, (right) Peak CO J= 2 — > 1 intensity along each line of sight for NGC 628. 
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of the telescope are the dominant sources of uncertainty. 
We estimate the uncertainty in the caUbration at ~ 20% 
(Section 13. 5p and find that applying masking (e.g., as in 
Figures [T7l - [30]) alters the flux by < 15% (slightly more 
in the cases of NGC 925, NGC 4214, and NGC 2841). 

In four low metallicity dwarf irregulars (Holmberg I, 
Holmberg II, IC 2574, and DDO 154) we do not detect an 
extended distribution of CO. Instead, Table[3]gives upper 
limits on their CO content. To derive these, we create 
an integrated spectrum for each cube and estimate the 
RMS noise at velocities offset from the H I line, regions 
expected to be signal-free. Combin ing this RMS noise 
with the H I hne width {W20 from IWalter et aP ^QQB) 
and area of the map yields the 5a upper limits on Fqo 
in Table [3l This approach takes into account the low- 
level correlated noise discussed in Section [23] 

We quote CO flux, Fco, in units of K km s~^ arcsec^ 
(Tmb)- This may be converted to Jy km via 

Fco [Jy km s"^] = 0.036 Fco [K km s~^ arcsec^] , 

(2) 

We report CO luminosity, icOj in units of K km s ^ pc^ 
(2mb)- This may be combined with a CO-to-H2 conver- 
sion factor to estimate the mass of molecular gas via 

Mh2 [Mq] = 5.5 X2 ^ Leo [K km s-^ pc^] . (3) 

Here X2 is the CO J ~ 1 ^ 0-to-H2 conversion fac- 
tor normalized to 2 x lO^*' cm^^ (K k m s~^)~^, approxi- 
mately th e Solar N eighborhood value (jStrong fc Matto^ 
[1996; Dani_eLalJ [2001). R21 is the CO J = 2 ^ 1 to 
J = 1 ^ ratio; 0.8 is a typical value that we find com- 
paring HERACLES to other surveys ( §5.4.2|) . Equation 
[3| includes a factor of 1.36 to account for helium. 

4.4. GMC Filling Fraction and Point Source Sensitivity 

In the Milky Way and nearby disk galaxies, CO emis- 
sion is observed to come mostly from GMCs. These 
clouds have approximately constant surface densities, 
Eh2 170 Mfn pc~^, and brightness temperatures > 4 K 
(jSolomon et al.lll987l l. We expect that these clouds sim- 
ilarly dominate the ISM across the disks of our sample. 

The contrast between our measured surface densities 
and those of Galactic GMCs may give an idea of the 
filling factor of GMCs inside a resolution element. To 
translate CO intensity, /cO: into molecular gas surface 
density, Eh2 , one can apply an analog of Equation [3l 

Eh2 [Mo pc^] = 5.5 X2 ^ Ico [K km s'^] . (4) 

The contours in Figures [iTj - [30] thus correspond to 
Sh2 « 5.5 (gray), 11 (first black), 22, 44, 88, 176 (first 
white) M0 pc~^, and so on, values that can be explained 
if ~ 3, 6, 13, 25, 50, and 100% of the area inside the beam 
is covered by Galactic GMCs. Most area where we see 
CO emission has relatively low filling factors, ~ 5-15%. 

The peak intensity maps allow a similar comparison. 
If ^ 4 K is the typical brightness for a Galactic cloud, 
Tpcak = 0.1 and 0.5 K — values that bracket most of 
the observed peak temperatures — correspond to filling 
factors of ~ 3 and 13%. 



The point mass sensitivity of our maps corresponds 
to a collection of several Galactic GMCs. For a source 
with full li ne width of ^ 10.4 km s~^ (typical for Galac- 
tic GMCs, [Solomon et"allll987f) . our typical RMS sen- 
sitivity translates to a 3(t point source sensitivity of 
4.8 X 10^ dio K km s~^ pc^, where diQ is the distance 
to the source divided by 10 Mpc. From Equation[31 this 
corresponds to an H2 mass of 2.7 x 10^ c?io Mq, which 
is near the high end found for Galactic GMCs (clouds 
with masses lOf-lOfM0 account for most of the Galac- 
tic H2, e.g.. [Blit3[1993( ). For reference, the most massive 
Milky Way GMC, the most massive cloud detected in 
M33, and the Orion molecular complex have masses of 
- 6 X 10^ Mq, ^ 7 X 10^ Mq and ^ 5 X 10^ Mp, re- 
spectively (Solomo n et all [19871 : [Engargiola et all 120031 : 
iWilson et al . 2005). 

5. COMPARISON WITH OTHER DATA 

Previous observations of our targets at wavelengths 
from radio to UV allow a range of comparisons. Here 
we use these data to calculate the i72-to-H I and i?2-to- 
stellar mass ratio, to test our assumption that H I and 
CO exhibit the same mean velocity, to check our obser- 
vations against previous CO measurements, and to esti- 
mate the CO J = 2 ^ 1/1^0 line ratio. More detailed 
comparisons to H I, IR and UV da ta may be found in 
[Lerov et al.l ([2008[ ) and lBigiel et"an ([2008D . 

5.1. Relative Masses of H2, H I, and Stars 

Columns 7, 8, and 9 in Table [3| list the ratios of H2 
mass to H I mass {M-ii2 / M-m) , H2 mass to stellar mass 
(Mh2/M*), and gas to stellar mass ((Mh2 + Mm) /M^) 
for ea ch galaxy. We adopt Mhi from [Walter et al.l 
([2008f ). scaling by 1.36 to a ccount for h e lium. The 
stellar masses are taken from [Lerov et al.[ ([2008[ ) , who 
estimate them from SINGS 3.6/xm ima ging assuming 
a Kr ou pa IMF (NGC 2903 comes from Ide Blok et al.l 
[2008D . Ide Blok et all ([2008D investigate possible varia- 
tions in the 3.6/Ltm mass-to-light ratio and find a factor 
of 2 plausible. 

The range of A/H2/Af* is 0.01-0.25, with most values 
between 0.03 and 0.1. A/h2/-^^hi spans from 0.02 to 1.13, 
with most targets in the range ~ 0.2-0.6 (in our sam- 
ple only NGC 4736 appears to have more H2 than H I) . 
Combined, these ratios yield a wide range of gas-richness, 
from {Aim + M}i2) /M^, ~ 0.2 in relatively early type spi- 
rals to > 3 in our late-type, low-mass nondetections. 

5.2. Scale Lengths At Other Wavelengths 

Figure [32] shows the size of our galaxies at other wave- 
lengths as a function of the exponential scale lengths that 
we measure from azimuthally avera ged CO profil e s, Icn - 
Filled circles show 0.2 r25, which [Young et al] ([1995] ) 
found to be a typical scale length of CO emission. We 
observe the same here, Ico = 0.2 r25 with a la scat- 
ter of 0.05. Stars show the exponential sca le length fit 
to m edian profiles of 3.6/im emission (from ITTerov et ahl 
[2008f ). a tracer of the distribution of old stellar mass. 
Diamonds plot the exponential scale length fit to pro- 
files of star formation surface density, estimated from a 
combinatio n of GALEX FVY emission and Spitzer 24/im 
data (from lLeroy et al]l2008l ). All three quantities scat- 
ter about equality, with Ico on average ~ 5-10% lower 
than the other two, not a significant difference. The 
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Fig. 18.— As Figure [17] but for NGC 925. 
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Fig. 20.— As Figure [TtI but for NGC 2903. 
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Fig. 21.— As Figure [TTl but for NGC 2976. 
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Fig. 23.— As Figure [TTl but for NGC 3198. 
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Fig. 24.— As Figure [TTl but for NGC 3351. 

plot underscores the weU-established close association be- 
tween CO emission and stellar light in disk galaxies (e.g., 
Young fc Scovilld 119911 : lYoung et aH ll995HRegan et all 
2001[ ) and the similar match between the distrib utions of 
CO and star formation fe.g.. iLerov et al.ll2008[ ). 

5.3. HERACLES CO and THINGS Hi Velocities 

The THINGS data allow us to test our assumption 
that CO and H I have the same mean velocity fi i3.2[l . To 
do so, we assemble a set of ^ 1200 spatially indepen- 
dent spectra, each with peak SNR greater than 5. For 
each spectrum, we measure the intensity-weighted mean 
velocity of CO J = 2 ^ 1 emission, {vco), and extract 
the intensity- weighted mean H I velocity, (whi), from the 
THINGS first moment maps (Walter et al. 2008). 

The gray histogram in Figure [33] shows the distribution 
of {vco) — {vm) for all 1200 spectra. The mean velocities 
of CO and H I appear closely matched for most spectra. 



A Gaussian fit to the distribution (thick dashed line) has 
Icr width 5.9 km and a center consistent with zero off- 
set (—0.05 km s~^). The scatter in (uco) — (vm) is much 
smaller than the 60-300 km s~^ width of our baseline fit- 
ting regions (i j3.2p . so the agreement is not a product of 
our reduction. The extended wings come almost entirely 
from a few inclined, massive spirals (NGC 2841, 2903, 
3521, 7331), galaxies that also sometimes show multi- 
ple components in a single H I spectrum. These outliers 
arc interesting, but here we emphasize the basic result 
that H I and CO gas show approximately the same mean 
velocity. 

5.4. Comparison to CO J = 1 — > Observations 

CO J = 1 ^ emission has been mapped for many 
galaxies in our sampl e. Here we compar e HERACLES to 
the FCRAO s urvev (I Young et al.|[l995l resolution 45"), 
BIMA SONG (jHelfer et al.l2003f) . and the Nobeyama CO 
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Fig. 25.— As Figure [TTl but for NGC 3521. 
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Fig. 26.— As Figure [TTl but for NGC 4214. 
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Fig. 27.— As Figure [l3 but for NGC 4736. 




ff 60^46" 

R.A. (JSQDO) 



HERACLES: The HERA CO-Line Extragalactic Survey 



21 



20 in BO D.Q n.i o.s o.s 

IntefjrBted CO Intensity [K km a~'] Peak Intensity [E] 




Fig. 28.— As Figure fTTI but for NGC 5055. 
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Fig. 30.— As Figure [TtI but for NGC 7331. 

Atlas of Nearby Spiral Galaxies (|Kuno "eTall [2007 1 res- 
olution 15"). We use two data sets drawn from BIMA 
SONG: on-the-fly maps made with the NRAO 12-m (res- 
olution 55") and the combined BIMA -I- 12-m data (typ- 
ical resolution ^ 6"). We refer to the former as the 
"NRAO 12-m" and the latter as BIMA SONG. 
Our method of comparison is the following: 

1. We convolve each plane of each HERACLES data 
cube with a series of Gaussian kernels to create 
data cubes with angular resolutions of 15, 45, and 
55", a ppropriate for comparison with iKuno et al.l 
()2007D . the FCRAO survey, and the NRAO 12-m 
data. We convolve BIMA SONG to 15" resolution 
in the same manner and compare it to HERACLES 
at this resolution. 

2. For each galaxy, we extract spectra from each data 
cube at a series of independent pointings. When 
comparing to the FCRAO survey, these are the 
FCRAO pointings, which are usually spaced by a 
full (45") beam width along the m ajor axis. When 
comparing with lKuno et al.l (|2007f ). BIMA SONG, 
and the NRAO 12-m, the pointings are on a grid 
that covers most of the galaxy with sampling points 
spaced by a full beam width — i.e., 15" and 55". 

3. We smooth these spectra in velocity so that they all 
have comparable velocity resolution, ^ 15 km s^^ 
(chosen to match the FCRAO survey). 

4. We discard all spectra where the peak SNR is less 
than 5. When comparing with the FCRAO sur- 
vey, w e consider only pointings where lYoung et al.l 
(|1995f ) report a peak temperature, mean velocity, 
and line width. 



The result is a series of reasonably high SNR spectra 
at matched positions with matching angular and velocity 
resolutions. For each spectrum, we measure four param- 
eters: the peak intensity of the line, Tpeak (in K); the 
intensity-weighted mean velocity of the line, Umcan (in 
km s""'^); the full width at half maximum of the line, 
fpwHM (in km s~^); and the integrated intensity of the 
line, Jco (in K km s~^). Because the FCRAO survey 
data are not electronically av ailable, we use the v alues of 
these parameters reported bv lYoung et al.l (|1995f ). 

5.4.1. CO Velocity and Line Width 

The line width and mean velocity are not expected to 
vary strongly between the CO J = 2 ^ 1 and J = 1 ^ 
transition. Therefore comparing our measurements of 
these quantities to those from earlier surveys allows a 
basic check on our data. Overall, this exercise confirms 
that we measure the same basic line shapes found by pre- 
vious surveys at matched positions and resolution. We 
show this in the top two panels of Figure IMI and the first 
two lines of Table |4l which give mean velocity offset and 
line width ratio (with equal weight to each galaxy). 

The top left panel of Figure [M] shows the distribu- 
tion of differences between the intensity-weighted mean 
CO J = 2 ^ 1 (HERACLES) velocity and the intensity- 
weighted mean CO J = 1 ^ velocity. Each comparison 
survey appears as a separate, normalized histogram. A 
thick dashed line shows a Gaussian f it to the ave r age o f 
the BIMA SONG, NRAO 12-m, and lKuno etlHI ([20071 ) 
histograms. This Gaussian is centered at —1.2 km s~^ 
and has Icr width 5.6 km s~^, identical with the uncer- 
tainties to what we found comparing CO J = 2 — s- 1 to 
HI. 

The top right panel of Figure IMI shows the FWHM line 
width, upwHM, from HERACLES measured as a function 
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Black circles and gray triangles show the intensity averaged over 10"-wide tilted rings. Gray triangles indicate where the measurement 
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TABLE 4 

Comparison of HERACLES and Other Surveys 



Quantity FCRAO Survey NRAO 12-m BIMA SONG Kuno et al. (2007) 

Aumcan (km s"!)^ -2.9 3.6 2.4 -0.4 

ratio of i'fwhm'' l-Ol 1.05 1.18 0.93 

ratio of Tp^ak 

0.76 0.87 0.67 0.81 

ratio of /co'' 0.71 0.87 0.84 0.74 



^ Values are mean diflterence (HERACLES - other survey) or ratio (HERA- 
CLES/other survey) giving equal weight to each galaxy. 
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Fig. 32. — Several sizes for our targets (y-axis) as a function of 
the exponential scale lengths fit to the azimuthally average pro- 
files (x-axis). Circles show 0.2 r25; stars show the 3.6/jm scale 
length, a tracer of the old stellar distribution; and diamonds show 
the scale length GALEX FUV -I- Spitzer 24/xm emission. Vertical 
lines connect points from the same galaxy. The dashed line shows 
equality. 



of upwHM from other surveys. Error bars show la un- 
certainty, estimated by repeatedly adding the measured 
noise to each spectrum and re-measuring wfwhm (with- 
out the FCRAO spectra, we cannot estimate uncertain- 
ties for these data). As with the mean velocity, there 
is general good agreement be tween HE RACLES, BIMA 
SONG, the NRAO 12-m, and lKuno et all (|2007h . There 
are systematic differences: BIMA SONG tends to show 
sUg htly lower line wid ths than the other data sets while 
the iKuno et al.l (|2007f) yield slightly higher line widths. 
These have magnitude ±20% and affect the intercompar- 
ison of CO J = 1 ^ data as much as the comparison 
between HERACLES and the other surveys. 

The outlier this comparison is the FCRAO survey. 
These data do not agree as well with our own as the other 
three data sets. Specifically, there are several pointings 
where the mean velocity and line width disagree strongly 
with our data. Most of these loca ti ons a lso overlap with 
BIMA SONG or the iKuno et al.l (|2007f ) survey and in 
these cases, these surveys also disagree with the FCRAO 
survey. We note that 1) the FCRAO survey is the only 
data set from which we do not ourselves measure the line 
parameters, so methodology may drive some of the dif- 
ference; 2) these are the only data which are not maps, so 
it is not easy to tell whether pointing offsets may affect 
the comparison. 
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Fig. 33. — Distribution of (vqq) — (dhi)i thi^ difference between 
the intensity-weighted mean CO velocity (from HERACLES) and 
intensity- weighted mean H I velocity (from THINGS), along lines 
of sight with good signal-to-noise. The best fit Gaussian (dashed 
line) has lo" width 5.9 km and center —0.1 km s~^. 



5.4.2. CO J = 2^1/J =1^0 Lme Ratio 

The ratio of CO J = 2 ^ 1 measured by HERACLES 
to CO J = 1 ^ emission measured by previous sur- 
veys reflects both the relative calibrations of the tele- 
scopes and the line ratio, R21 = Ijtwo/Ii — > Q, which 
depends on the optical depth and excitation tempera- 
ture of the gas. Generally speaking, R21 > 1 indicates 
warm, optically thin gas. Optically thick gas produces 
i?2i ~ 0.5-1.0 with the value reflecting the temperature 
of the gas, though in low density regions the levels may 
not be populated accordi ng to local thermod ynamic equi- 
librium (LTE; see, e.g., l Eckart et al.lll990l . for a more 
thorough discussion). 

In the bottom left panel of Figure [M] we plot the 
peak temperatures measured by HERACLES (y-axis) as 
a function of the peak temperature measured from other 
surveys at matched position and resolution (a;-axis). We 
use the peak temperature because it avoids any concern 
related to the line profile or range of integration. We 
do not plot error bars, but recall that our condition to 
include a spectrum in this analysis was that Tpcak have 
SNR > 5 (meaning the statistical uncertainty on the ra- 
tio is always la < 30%). Table [5] compiles the average 
peak temperature ratio for each galaxy and survey and 
notes several measurements from the literature. 
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TABLE 5 

Average Tpcak Ratio by Galaxy and Survey^ 



Galaxy 


FCRAO Survey 


NRAO 12-m 


BIMA SONG 


Kuno et al. ("2007") 


Other*" 


NGC 628 


0.93 


1.01 






0.44 


NGC 2841 


0.72 










NGC 2903 


0.91 


0.90 


0.68 


0.82 


0.6-0.9 


NGC 3184 


0.59 






0.60 




NGC 3351 


0.50 


0.62 


0.59 


0.95 


1.65 


NGC 3521 


1.14 


0.81 


0.59 


0.86 




NGC 4736 


0.70 


1.02 


0.83 


1.12 




NGC 5055 


0.71 


0.77 


0.64 


0.79 




NGC 6946 


0.85 


1.03 


0.77 


0.92 


0.8 


NGC 7331 


0.82 


0.79 


0.60 




0.5 



^ Mean ratio of peak CO J — 2^1 temperature to peak CO J = 1^0 
temperature.'' Integrated intensity {not peak temperature) ra tios drawn from th e 
literature. R eferences: NGC 628 NGC 3351, NGC 7331 — iBraine et all (Il993l ): 
NGC 2903 — iJackson eraP l|1991f ): NGC 6946 — ICrosthwaite fc Turned (|2007f i 



The ratios in Table [5] span 0.48-1.06, with most val- 
ues between 0.6 and 1.0. These values agree with pre- 
vious observations of the Milky Way and other galax- 
ies. CO emission in the Milky Way comes mostly from 
optically thick regions and typical line ratios are ~ 0.6 
- 0. 8 in the disk and ~ 1 towards the Galactic cen- 
ter (jSakamoto et al.l 119951: lOka et al.l [19 96. 199 81). In 
the ce ntral regions of other galaxies. iBraine fc Combed 



1992 ) found an average line ratio ^ 0.9. ICasoli et al.l 



199lh compile observations on a number of galaxies and 
report typical values of ~ 1 in galaxy nuclei and 0.5 - 
0.7 in galaxy disks. 

If the gas is optically thick, a ratio of 0.6 corresponds 
to an excitation temperature of ^ 5 K, 0.8 to ~ 10 K, 
and 0.9 to ^ 21 K. Higher excitation temperatures yield 
a line ratio ~ 1. The average of all the values in Ta- 
ble E] is 0.81, which is consistent with optically thick 
gas with an excitation temperature ~ 10 K. One should 
treat this temperature with caution: we do not constrain 
the optical depth (e.g., using one of the isotopes of CO) 
and we have no verification that our sources are in LTE, 
so sub-thermally excited, low-density envelopes — which 
are observed in Gal actic GMCs — may low er the overall 
intensity ratio fe.g.. ISakamoto et al.lll994f) . 

There is noticeable scatter among the peak tempera- 
tures ratios derived from different data sets for the same 
galaxy. Some scatter may arise from differences in reso- 
lutions or area considered (the set of spectra with peak 
SNR > 5 varies with comparison survey). If we ignore 
these explanations. Table [Hj suggests that there are sys- 
tematic calibration differences among the various surveys 
at the ~ 10% level and that the calibration of a given 
galaxy in a given survey is uncertain by ~ 15%. These 
uncertainties seem reasonable for millimeter line observa- 
tions and agree with the scatte r in comparisons w ith the 
FCRAO surve y carried oui^ bylHelfer et al.l (j2003l . their 
Figure 51) and lKuno et al.l (|2007l . their Figure 1). 

Because R21 is observed to vary between the center 
and the disk in both the Milky Way and other galaxies, 
the bottom right panel of Figure |M| shows the ratio of 
Tpeak in the two transitions as a function of galactocen- 
tric radius. This plot s hows only data from HERACLES 
and iKuno et al.l (j2007D . These data most closely match 
our own in observing strategy, instrumentation, and res- 
olution making them the best option for this comparison. 



Black points connected by lines show the median ratio 
in bins 0.05 r25 with the first bin centered at 0. In the 
centers of galaxies , we observe the same trend noted by 
ICasoh et al] ()199lD and found in the Milky Way: the ra- 
tio is high 1.3) in the center of the galaxy and then 
drops rapidly to ~ 0.8, remaining almost constant at this 
value out to the edge of our comparison at ^ 0.45 r25. 
This constancy must be interpreted with care. By se- 
lecting only high significance spectra, we bias ourselves 
to bright regions. This may have the effect of creating 
a homogeneous data set while omitting lower intensity 
emission that accounts for a significant fraction of the 
emission at large radii. 

6. SUMMARY 

This paper presents the HERA CO-Line Extragalactic 
Survey (HERACLES), a new survey of CO J = 2 ^ 1 
emission from nearby galaxies obtained using the IRAM 
30m. Because on-the-fly mapping mode with the multi- 
pixel HERA allows us to efficiently map large areas with 
good sensitivity and adequate resolution, HERACLES 
targets a wider area than previous surveys. Most maps 
extend to the edge of the optical disk, defined by the 
isophotal radius r25. 

Our sample overlaps THINGS and SINGS, allowing 
easy comparison to data from radio to UV wavelengths. 
One application of this multiwavelength data allows us 
to automate our reduction. We use the mean H I velocity 
from THINGS as a prior to define baseline fitting regions 
and are thus able to reduce 10^ spectra in a robust, auto- 
mated manner. We verify the assumption that CO and 
H I exhibit the same mean velocity by direct comparison 
of HERACLES and THINGS mean velocities, finding ht- 
tle or no systematic difference. We also confirm that our 
automated approach matches a by-hand reduction well. 

We clearly detect 14 galaxies and place upper limits 
on CO jtwo emission from 4 low metallicity, irregular 
dwarf galaxies. For a Galactic CO-to-H2 conversion fac- 
tor, the implied H2 masses are mostly 0.2-0.6 times the 
H I mass and 0.03-0.1 times the stellar mass. More de- 
tailed analysis of the relationship between H2, HI, stars, 
and star formation in our sample is presented elsewhere 
(Bigiel et al. 2008; Lerov et al.l[2008f ). 

We illustrate the distribution of CO in the detected 
galaxies via channel maps, azimuthally averaged profiles, 
and intensity maps. Where we find emission, the bright- 
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intensity weighted mean velocity measured by HERACLES a nd four surveys of C O 1 — > emission. A thick dashed line shows a Gaussian 
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indicates an individual pair of high significance spectra. The black points and lines show the median and scatter for the data in bins 
0.05 r25 wide with the first bin centered at 0. 



ness is usually consistent with ~ 5-15% of the area inside 
the beam being covered by Galactic GMCs, though in a 
few regions the implied surface density averaged over our 
^ 500 pc beam exceeds that of a Galactic GMC. The line 
width and mean velocity that we derive agree reasonably 
with previous (CO J = 1 — s- 0) surveys; systematic dif- 
ferences among surveys do exist, with magnitude '--^ 10- 
20%, but HERACLES never appears to be an outlier. 

The ratio of CO J = 2 ^ 1 intensity to CO J = 1 ^ 
intensity for high significance spectra lies mostly in the 
range 0.6 - 1.0 with an average value of 0.8, comparable 
to that found in the inner Milky Way. This could be 
produced by optically thick gas with an excitation tem- 
perature 10 K, though this temperature should not be 



over-interpreted without constraints on the optical depth 
or applicability of LTE. This line ratio is higher (~ 1.3) 
in the centers of galaxies and then roughly constant as a 
function of radius, though we caution that strong selec- 
tion effects may be at work. 

Our detections include high significance emission from 
the outer part of the star-forming disk in several galax- 
ies. We also find bright CO emission associated with an 
H I filament outside the main disk of NGC 4736 (radius 
1.33 r25). The azimuthally averaged intensity usually 
declines smoothly even beyond these detections and we 
are sensitive to only the most massive individual GMCs. 
Therefore we would expect to find many more such re- 
gions with improved sensitivity. 
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Although the channel and integrated intensity maps 
show a variety of morphologies, the azimuthally averaged 
profiles are well-described by exponential declines. The 
best-fit scale lengths range from 0.8 - 3.2 kpc and cor- 
relate closely with optical radius, near-IR (stellar mass) 
scale length, and UV-I-IR (star formation) scale length. 
These results are in good agreement with previous find- 
ings that CO emission closely follows both the stellar 
light and distribution of star formation on large scales. 
The exponential decline in CO brightness is a combina- 
tion of decreasing maximum brightness of CO emission 
and a decreased filling fraction of bright CO emission at 
large radii. 
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